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Abstract 

The Fresnel equation governing the propagation of electromagnetic waves for 
the most general linear constitutive law is derived. The wave normals are 
found to lie, in general, on a fourth order surface. When the constitutive 
coefficients satisfy the so-called reciprocity or closure relation, one can define 
a duality operator on the space of the two-forms. We prove that the closure 
relation is a sufficient condition for the reduction of the fourth order surface 
to the familiar second order light cone structure. We finally study whether 
this condition is also necessary. 
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I. INTRODUCTION 



Electromagnetic wave represents perhaps the most important classical device with the 
help of which one can carry out physical measurements and transmit information. Intrinsic 
properties and motion of material media, as well as the geometrical structure of space- 
time, can affect the propagation of electromagnetic waves. In the most general setting 
|T|J^, electromagnetic phenomena are described by the pair of two- forms H,F (called the 
electromagnetic excitation and the field strength, respectively) which satisfy the Maxwell 
equations dH = J,dF = 0, together with the constitutive law H = H{F). The latter 
relation contains a crucial information about the underlying physical continuum (i.e., mate- 
rial medium and/or spacetime). Mathematically, this constitutive law arises either from a 
suitable phenomenological theory of a medium or from the electromagnetic field Lagrangian. 

In general, the constitutive law establishes a nonlinear (or even nonlocal) relation be- 
tween the electromagnetic excitation and the field strength. The function (or functional) 
H{F) may depend on the polarization and magnetization properties of matter, and/or on the 
spacetime geometry, i.e. metric, curvature, torsion and nonmetricity. Previously, propaga- 
tion of electromagnetic waves was analysed for a variety of constitutive laws: for nonlinear 
models in Minkowski and Riemannian spacetimes [^, for electrodynamics in a Riemann- 
Cartan manifold [Q, and also for certain nonminimal and higher derivative gravity models 
1^. Numerous authors |^ discussed the electromagnetic waves in the Einstein-Maxwell the- 
ory. The main aim of this paper is to investigate the wave propagation in the Maxwell 
electrodynamics with the most general linear constitutive law. We derive the generalized 
Fresnel equation which determines the wave normals directly from the constitutive coeffi- 
cients. This result is of interest, e.g., for various applications in crystaloptics and related 
domains. 

Another motivation for the present work came from the study of a deep relationship be- 
tween the duality operators defined on two-forms and the conformal classes of spacetime met- 
rics in four dimensions. Within the classical Maxwell electrodynamics, Toupin, Schonberg 
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and others ||^ have noticed that the constitutive coefficients define a duahty operator, pro- 
vided a certain reciprocity or closure condition is fuffilled, and gave first demonstrations 
of the existence of the corresponding conformal metric structure. Later these observations 
were rediscovered and developed in mathematics [H and in the gravity theory [ffO]. Recently 



the complete explicit solution of the closure relation has been given []TT|, and it was conjec- 
tured that the reciprocity condition is a necessary and sufficient condition for the standard 



null-cone structure for the light propagation (see also independent arguments in [|T3I)- Here 
we give a partial answer to this question. 



II. ELECTRODYNAMICS WITH LINEAR CONSTITUTIVE LAW 

Let us consider the Maxwell equations in vacuum, 

dH = 0, dF = 0, (2.1) 

i.e. we assume that the electric current three-form J vanishes in the spacetime region under 
consideration. Given the local coordinates x*, i = 0,1, 2, 3, we can decompose the exterior 
forms as 

H = ^ Hij dx' A dx\ F=^ Fij dx' A dx^ . (2.2) 

Following [|rT],|T3|, we write the linear constitutive law in terms of the electromagnetic exci- 
tation and field strength tensors as 

= ^ e,,M x''"^" F^r^, ^, J, . . . = 0, 1, 2, 3. (2.3) 

Here eijki is the Levi-Civita symbol and x*-'^'(x) an even tensor density of weight +1 (called 
the constitutive tensor density) which can be decomposed according to 

^ijki ^ j^^^ ° i,ki ^ ^^^^ ^ijki ^ ^.^j^ ° [i,ki] ^ Q _ (2.4) 

Here f{x) is a dimensionfull scalar function such that X ^-^ is dimensionless. The pseudo- 
scalar constitutive function a{x) can be identified (on the kinematic level) as an Abelian 
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axion field, whereas /(x) can be interpreted as a dilaton scalar field. Note that X*-' has the 
same algebraic symmetries and therefore the same number of 20 independent components 
ClS db Riemannian curvature tensor: 



^ ijkl _ _ ^ jikl _ _ ^ ijlk _ ^ klij ^ [ijkl] _ g 



(2.5) 



This follows from the existence and the structure of the Lagrangian for the linear electro- 
dynamics Viin = — \H/\F, see P,p!3||. It is convenient to adopt a more compact (essentially 
bivector) notation by defining the 3- (co) vector quantities 



-noi 



02 



and 



F31 



Ea-.-- 



20 



y -^30 j 



(2.6) 



for the electric and magnetic excitations, and for the magnetic and electric field strengths, 
respectively. The Latin indices label now a, 6, c, . . . = 1, 2, 3. The constitutive tensor is then 
naturally parametrized by a triplet of 3 x 3 matrices, T^^' = {^«^S.,,C^}, so that the 
constitutive law ( p. 41 ) is finally recasted into 



( nb n \ ^ 



I 



+ a[x) 



(2.7) 



Here the 3x3 matrices satisfy A 



ah 



ah 



Bba, and C^a = 0, thereby providing the 



algebraic properties (|2.5|). 



III. WAVE PROPAGATION: FRESNEL EQUATION 

In the theory of partial differential equations, the propagation of waves is described by 
Hadamard discontinuities of solutions across a characteristic (wave front) hypersurface 5" 
0. One can locally define S by the equation $(x*) = const. The Hadamard discontinuity 
of any function jF(x) across the hypersurface S is determined as the difference [JF] (x) := 
), where x± := lim (x ± e) are points on the opposite sides of S 3 x. An 

e— >0 
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ordinary electromagnetic wave is a solution of the Maxwell equations (|2.1|) for which the 
derivatives of H and F have regular discontinuities across the wave front hypersurface S. 
In terms of the (co)vector components, we have on the characteristic hypersurface S: 



= 0, [d^V^] = d'' q,, [Ha] = 0, [d^Ha] = K Qi 

[B^] = 0, [d,B^] = b'' g„ [Ea] = 0, [diEa] = g„ 



(3.1) 
(3.2) 



where d°',ha,b°', Cq describe discontinuities of the corresponding quantities across S, and the 
wave-covector normal to the front is given by 



Qi := di^ 



(3.3) 



Equations (|3.1| )- (|3^ ) represent the Hadamard geometrical compatibility conditions. Substi- 



tuting (13) into dOI), and using (|J) and (|]I])-(g]2D, we find 



go d"" - e""'' qt K = 0, go + e"^"^ = 0, 
qad'' = 0, qab'' = 0, 



(3.4) 
(3.5) 



where e"'"^ is the 3-dimensional Levi-Civita symbol. In this system only the 6 equations 
( ^.41 ) are independent. Assuming that go 7^ 0, one finds that the equations ( |3.5D are trivially 
satisfied if one substitutes (|3.4| ) into them. [Note that the characteristics with go = do 
not have intrinsic meaning for the evolution equations, since they obviously depend on the 
arbitrary choice of coordinates]. 

Differentiating ( p.7| ) and using the compatibility conditions (|3.1| ) -(|3.2| ) , we find addition- 
ally 6 algebraic equations: 



= /(^) 

V 



C a Bab 



l-e ^ 



(3.6) 



Note that the constitutive coefficients and their first derivatives are assumed to be continuous 
across S. 

We can now substitute d"" and ha from (|3.6| ) into the first equation ( p.4|) , which gives 



f{x)qo {-A''' e, + C\ b') + a{x)qo = fix^'^qb {-C\ + 3,^ h") - a[xY''-q^ e,. (3.7) 

The terms proportional to the axion field a{x) drop out completely due to ( |3.4| ), and then 
one can also remove the common dilaton factor f{x) on both sides of the equation. [We 
assume f{x) 7^ 0, since otherwise there is no hyperbolic evolution system]. It remains finally 
to substitute 6" in terms of e;, from the second equation ( p.4|) , and after some rearrangements 
one finds: 



2 jiab I \ no- cdb , /^b cdaA , aec bfd f2 \ n 

q^A + qoqd C +C ce + qeqj e e^Bcd)eb = 



This homogeneous algebraic equation has a nontrivial solution when 



W := det 



„2 Aab I „ „ ^cdb , nb ^cda , „ „ nechfd to 

QqA + qoqd t-cC +(-^ + Qeqj e e t^cd 



0. 



(3.8) 



(3.9) 



This is a Fresnel equation which is central in the wave propagation analysis. It determines 
the geometry of the wave normals in terms of the constitutive coefficients A, B, C. A direct 
calculation yields the general result: 

W = ql (q^M + qlqa + qlq^qb M"^ + q^qaqtqc M"''^ + qaqbqcqd M''""') = 0, (3.10) 
where we have denoted 



M := det A 



■.= 2e,,dA^'C%A^'', 



M"' := BcdiA""' A"" - A' A'") - A"" C\ C\ + 4 A' C\ C'c - C'd 



abc 



2e 



cde 



BdfiA'^'' - A^^ C\) + C\ C'f C^, 



M' 



abed 



.— ^cef^dgh ^ 



fh 
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(3.11) 
(3.12) 
(3.13) 
(3.14) 



Note that only the completely symmetric parts M*^°i'"'^p\ p = 2,3,4, contribute to the 
Fresnel equation. Since go 7^ 0, one can delete the first factor in (|3.10|) , and thus we finally 
find that the wave covector lies, in general, on a 4th order surface. This is different from 
the light cone (i.e., 2nd order) structure which arises only in a particular case. In the next 
section we demonstrate that the latter corresponds to the closure condition. Earlier, the 



relation between the fourth- and the second-order wave geometry was studied by Tamm 
for a special case of the linear constitutive law. 



IV. THE CLOSURE RELATION AS A SUFFICIENT CONDITION 



The linear constitutive law defines a duality operator when the constitutive coefficients 
satisfy the 'reciprocity' or 'closure' relation 
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(4.1) 



or in terms of the 3x3 matrices: 



(4.2) 



The general solution of the closure condition (^3 

1 



reads 1 1 



A 



detS 

C°'b = B"''^ edbc k 



{k'^B''^ -k" k^)-B''\ 
1 



(4.3) 
(4.4) 



deti3 

Here k"" is an arbitrary 3-vector, k^ := Babk"', k"^ := Babk°'k'', and B"'^ denotes the inverse 
matrix to Bab- 



Starting from (|4.3| )- (|4.4| ), the direct calculation yields: 



M = 



detB 
1 



detS 



4fcM 1 - 



det-B; 



detB. 



M 



ab 



1 / h'2 

Ak^'k^ + 2i3'^^ 1 - 



deti3 



deti3 



(4.5) 

(4.6) 

(4.7) 
(4.8) 
(4.9) 



Substituting all this into the general Fresnel equation ( p.lQ ), we find 

k^ \ 2qoiqak^) 



1 



deti3| 

2 



det-B 



det;B| 



det B\ {qaqbB 



'ab\ 



(4.10) 



Here a = sign{det B), and g^^ is the (inverse) 4-dimensional metric which arises from the 
duality operator and the closure relation |1Tll , [T3| : 



^°'^ = -^^, (4.12) 



g"' = -^\deiB\B''\ (4.13) 

This metric (^jj (defined up to a conformal factor) always lias the Lorentzian signature, 
although it is not necessarily interpretable as a spacetime metric (this is a so called optical 
metric, in general; see, e.g., [1^). As shown in [0, the constitutive tensor density ( p.4| ) can 
be rewritten in terms of this metric as 

ijkl _ f /„\ ^ /— : ( Jk jl _ „ik il\ _|_ Jjkl 



x''"' = fix) [g^'g^' - g''g'') + aix) e'^'^ . (4.14) 

Thus we indeed recover the null cone g^g* = qiQj g"^^ = structure for the propagation of 
electromagnetic waves from our general analysis: provided the constitutive matrices satisfy 
the closure relation ([4.1|) -(^I^), the quartic surface (|3.10|) degenerates to the null cone for 



the induced metric gij. 



It is worthwhile to note that the Fresnel equation (|3.1CI| ) can be rewritten in an explicitly 
covariant form 



G'^'^qirnkqi = 0, tj,... = 0, 1, 2, 3, (4.15) 

where the fourth order totally symmetric tensor density G*-^*^' is constructed as the cubic 
polynomial of the components of the constitutive tensor: 



^^,H _ ^ ^jMk ^i)stu (4.16) 



[Here the total symmetrization is extended only over the four indices i,j,k,l with all the 



summation indices excluded]. Tamm JTgI has introduced a similar 'fourth-order metric' for 



the particular case of the linear constitutive law. 
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V. THE CLOSURE RELATION AS A NECESSARY CONDITION 



It was conjectured [11,13| that the closure relation is not only sufficient, but also a 
necessary condition for the reduction of the quartic geometry ( ^.lUI ) to the null cone. The 
complete proof of this conjecture requires a rather lengthy algebra and will be considered 
elsewhere. Here we demonstrate the validity of the necessary condition in a particular case 
when the matrix C = 0. 



Putting Cb = 0, we find from (|3TT| )- (|3l^) that = and M'^'"^ = 0, whereas 



M"^ = BcdiA""^ A"^ - A"' A'"^), (5.1) 



Consequently, (|3.1CI| ) reduces to 

W = q^(detAq^ + ql-f + detBa(3), (5.3) 

where a := A°'^qaqb, P := B°-^qaqb, and 7 := M°'''qaqb- Assuming that the last equation 
describes a null cone, one concludes that the roots for q^ should coincide and thus necessarily 

7^ = 4det^deti3a/3. (5.4) 

Let us write (det^det;B) = s \ det^det;B|, with s = sz(7r;,(det ^ det ^B). Then ( p^) yields 

2x/|det^detS| - = s X, 2./| det ^det i3| - = (5.5) 

V V 7 A 

where A is an arbitrary scalar factor. Recalling the definitions of 7, we then find 

A^^ = sX^B''\ (5.6) 
Consequently, M = det A = s A^/ det B and M"'' = 2X^B°''^, and therefore one verifies that 

^ = ^ i>^yo + s Qaqb B'^" det bY . (5.7) 



deti3 

We immediately see that for s = — 1 the quadratic form in ( |5.7| ) can have either the (H ) 

signature, or (+ + H — ). Similarly, for s = 1 the signature is either (+ + ++), or (+ H ). 



Therefore, the Fresnel equation describes a correct hght cone (hyperbohc) structure only in 
the case s = —1. 

Finally, one can verify that the above solutions satisfy 



1 



ijmn^pqrs 



which for s = — 1 reproduces the closure relation (|4.1| ) after a trivial rescaling of the con- 
stitutive tensor density (and subsequently absorbing the factor A into the 'dilaton' field 
/)• 

VI. CONCLUSIONS 

In this paper we have derived, extending the earlier results (see e.g., P, [14| , [T6| ) , the 
Fresnel equation governing the propagation of electromagnetic waves for the most general 
linear constitutive law. The wave covector lies, in general, on a fourth order surface. Such 
generic fourth order structure is not affected by the axion- and dilaton-like parts of the 
constitutive tensor. Note however that the linear constitutive law H = a{x)F does not lead 
to hyperbolic evolution equations, and hence necessarily /(x) 7^ 0. 

We have proved that the closure relation (|4.1j ) is a sufficient condition for the reduction of 
the fourth order surface to the familiar second order light cone structure. The corresponding 
family of conformally related metrics g coincides with that derived in see also ||13|| . 



This result may be considered as an alternative (as compared to Urbantke's scheme P,|TO|) 
derivation of the Lorentzian metric g from a duality operator. In terms of the Lagrangian, the 
closure relation is equivalent to the statement that V^un = — | [f (x) F A * F + a{x) F A F], 
where the Hodge operator * is defined by the metric g. 

For the special case C^b = we have proved that the requirement of reduction of the 
fourth order Fresnel structure to a second order one implies a relation between the consti- 
tutive coefficients which is slightly weaker than the closure relation ( |4.1| ), in that it allows 
for an arbitrary scalar factor. The latter though can be removed by the redefinition of the 
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dilaton field f{x). Also the signature of the resulting quadratic form is not fixed, so that 
one has to impose hyperbolicity as a separate condition. 

It is worthwhile to note that the results obtained can be directly applied to the refinement 
and generalization of the previous analyses of the observational tests of the equivalence 
principle. See, for instance, fTH], where some particular cases of the Fresnel equation have 



been studied in this context. 
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